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Abstract 



Quasi-degenerate neutrino mass models (QDN) which can explain the current 
data on neutrino masses and mixings, are studied. In the first part, we 
study the effect of CP-phases on QDN mass matrix (m LL ) obeying \i — r 
symmetry in normal hierarchical (QD-NH) and inverted hierarchical (QD- 
IH) patterns. The numerical predictions are consistent with observed data on 
(i) solar mixing angle (#12) which lies below tri-bimaximal (TBM) value, (ii) 
absolute neutrino masses consistent with 0^(3(3 decay mass parameter (m ee ) 
and (iii) cosmo logical upper bound Yl\ m i- m LL is parameterized using only 
two unknown parameters (e, rj) within fi — r symmetry. The second part deals 
with the estimation of observed baryon asymmetry of the universe (BAU) 
where we consider the Majorana CP violating phases (a, (3) and the Dirac 
neutrino mass matrix (miR)- mm is taken as either the charged lepton or the 
up quark mass matrix, a, (3 is derived from the heavy right-handed Majorana 
mass matrix Mrr. Mrr is generated from mn and mm through inversion 
of Type-I seesaw formula. The predictions for BAU are nearly consistent 
with observations for flavoured thermal leptogenesis scenario for Type-IA in 
both QD-NH and QD-IH models. We also observe some enhancement effects 
in flavour leptogenesis compared to non-flavour leptogenesis by a magnitude 
of order one. In non-thermal leptogenesis QD-NH Type-IA is the only model 
consistent with observed data on baryon asymmetry. QD-NH model appears 
to be more favourable than those of QD-IH. The predicted inflaton mass 
needed to produce the BAU is found to be ~ 10 10 GeV corresponding to 
the reheating temperature Tr = 10 6 GeV. The present analysis shows that 
the three absolute neutrino masses may exhibit quasi-degenerate pattern in 
nature. 
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1 Introduction 



The presently available tightest cosmological upper bound of the sum of 
three absolute neutrino masses, has come down to the lowest value, J2i m i < 
0.28 eV [1], and a larger value of neutrino mass m 8 > 0.1 eV in Quasi- 
Degenerate Neutrino (QDN) mass models, has therefore to be ruled out. 
Furthermore, the upper bound value of neutrino mass parameter \m ee \ < 0.27 
eV appeared in the neutrinoless double beta decay (0v(3(3) experiments [2], 
also disfavours larger values of neutrino mass eigenvalues with same CP- 
parity. Investigations on QDN models in normal hierarchical (QD-NH) and 
inverted hierarchical (QD-IH) patterns of the three absolute neutrino masses, 
require a detailed numerical analysis to check whether such QDN models 
can really accommodate lower values of absolute neutrino masses \rrii\ < 
0.09 eV which are consistent with the above cosmological bound [1]. In the 
next step, the QDN models are again applied for the prediction of baryon 
asymmetry (r] B ) of the universe [3]. In order to estimate the observed baryon 
asymmmetry t)b = (6.51q;5)10 -10 [3] from a given neutrino mass model, one 
usually starts with a suitable texture of light Majorana neutrino mass matrix 
{jnLL) and then relates it with the heavy Majorana neutrinos matrix (Mrr) 
and the Dirac neutrino mass matrix (m^) through the inversion of Type-I 
seesaw mechanism in an elegant way. Since the stucture of Yukawa matrix 
for Dirac neutrino is not known, we consider the texture of Dirac neutrino 
mass matrix tulr as either the charged lepton mass matrix or up quark mass 
matrix, as allowed by SO(10) GUT models for phenomenological analysis. 

In many of the theoretical calculations on leptogenesis, the complex CP 
violating phases are usually derived from the Majorana phases appearing in 
PMNS leptonic mixing matrix Upmns which diagonalizes tull. Hence in 
such approach tjill is no longer taken as the starting point. However, in 
the present work, we consider a different route for the origin of complex CP 
violating phases which are derived from Mrr in the estimation of baryon 
asymmetry of the universe. This theoretical possibility is the main part of 
the present investigations. 

We first introduce a general classification of QDN models based on their 
CP-parity patterns in the three absolute neutrino masses, and then param- 
eterize the mass matrices using only two unknown parameters (e,rj), which 
can reproduce correct predictions on neutrino oscillation mass parameters 
and mixing angles, consistent with the latest observational data. Though 
such parameterization is intuitive, it is quite realistic for phenomenological 
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analysis. The estimations of baryon asymmetry of the universe in the light 
of thermal and non-thermal leptogenesis, may thus serve as an additional 
criteria to discriminate the correct pattern of neutrino mass models and also 
to shed light on the structure of unknown Dirac neutrino mass matrix. 

The paper is organised as follows. In section 2 we parameterize the neu- 
trino mass matrix. In section 3, numerical analysis and predictions on baryon 
asymmetry are outlined. Finally in section 4 we conclude with a summary. 
In the Appendix we briefly summarize the formalism for estimating the lep- 
ton asymmetry in thermal leptogenesis through out-of-equilibrium decay of 
the heavy right-handed Majorana neutrinos. This is followed by flavoured 
thermal leptogenesis and non-thermal leptogenesis. 



2 Parameterization and computations 

2.1 Parameterizations of neutrino mass matrix 

A general \i-t symmetric neutrino mass matrix [4,5] with its four unknown 
independent matrix elements, requires at least four independent equations 
for realistic numerical solution, 



m LL = 



I mn rn 12 m 12 \ 

m 12 m 22 m 23 
\ m 12 m 23 m 22 J 



(1) 



The three mass eigenvalues to, and solar mixing angle 9 12 , are given by 
mi = mn — V% tan#i 2 mi 2 , 
m 2 = mn + V2cot9i 2 mi 2 , 
m 3 = m 22 -m 23 . 

4. on 2y/2mi 2 

tan 26*12 = • (2) 

mn - m 22 - m 23 

The observed mass-squared differences are calculated as 



Am 2 



12 



2 

32 



m 3 — m 



(3) 



In the basis where charged lepton mass matrix is diagonal, we have the 
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leptonic mixing matrix, Upmns — U, where 



UpMNS — 



( cos 6 12 — sin 6 \2 \ 

sin 0i2 cos 0i2 1 

V2 V2 y/2 

sin 612 cos 6i2 1 

V \/2 V2 s/2 J 



(4) 



The neutrino mass parameter m ee in 0uf3f3 decay and the sum of the absolute 
neutrino masses in WMAP cosmological bound J2i m i, are given respectively 
by 

m ee = \miU^ + m 2 Ul 2 + m 3 U^ 3 \,m cosmos = m 1 + m 2 + m 3 . (5) 

A general classification for three-fold quasi-degenerate neutrino mass mod- 
els [5] with respect to Majorana CP-phases in their three mass eigenval- 
ues, is adopted here. Diagonalization of left-handed Majorana neutrino 
mass matrix hill in eq.(l) is given by m LL = UDU T , where U is the 
diagonalising matrix in eq.(4) and Diag=D(mi, m 2 e ia , m 3 e tl3 ) is the diago- 
nal matrix with two unknown Majorana phases (cx,j3). In the basis where 
charged lepton mass matrix is diagonal, the leptonic mixing matrix is given by 
U = Upmns [6]- We then adopt the following classification according to their 
CP-parity patterns in the mass eigenvalues rrii namely Type IA: (+— 1-) for 
D=Diag(mi, — m 2 ,m 3 ); Type IB: (+++) for D=Diag(mi, m 2 , m 3 ) and Type- 
IC: for (+- 1— ) for D=Diag(mi, m 2 , —rn 3 ) respectively. We now introduce the 
following parameterization for \x-r symmetric neutrino mass matrices ttill 
which satisfy the above classifications [4,5] and present a detailed numerical 
analysis. 



2.2 Numerical analysis and predictions 

For detailed numerical analysis we first choose the light Majorana neutrino 
mass matrix itlll presented in section 2.1. These mass matrices which obey 
ji — t symmetry [4], have the ability to deviate the solar mixing angle from 
that of TBM [7]. Next we estimate the numerical values of the three ab- 
solute neutrino masses. As discussed before, we need to introduce the neu- 
trino mass scale m 3 as input parameter in addition to the observed data 
[8] on solar and atmospheric neutrino mass-squared differences (Am 21 and 
| A TO32I). The three input parameters are fixed at the stage of predictions 
of neutrino mass parameters and mixing angles given in Tables 2 and 3. 
These results are consistent with the recent data on neutrino oscillation pa- 
rameters [8]. For numerical analysis we use the best-fit values of the global 
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input 


calculated 


QD-NH 


QD-IH 


m 3 


P 


mi m 2 


nil m 2 


0.40 
0.10 
0.08 


0.015 

0.24 

0.375 


0.39689 0.39699 
0.08674 0.08718 
0.06264 0.06326 


0.40289 0.40299 
0.11104 0.11136 
0.09340 0.09381 



Table 1: The absolute neutrino masses in eV, estimated from oscillation data, 
using calculated ip = 0.031667 as defined in the text. 



neutrino oscillation observational data [8] on solar and atmospheric neu- 
trino mass-squared differences. Amf 2 = (m| — m\) = 7.60 x !0~ 5 eV 2 and 
|Ara§ 2 | = \m\ — ml\ = 2 AO x 10~ 3 eV 2 , and define the following parameters 

p = ^ A ™H and ip = j^TT' wnere m 3 i s the input quantity allowed by the 

latest cosmological bound. For QDN in normal hierachy (QD-NH) pattern, 
the other two mass eigenvalues are estimated from, 

m 2 = m 3 ^l - p; m x = m 3 yjl - p(l + ip) (6) 

and for QDN in inverted hierarchy(lH-QD) the mass eigenvalues are ex- 
tracted from, 



m 2 = m 3 ^Jl + p;m 1 = m 3 yl + p(l - tp). (7) 

For suitable input value of m 3 , one can estimate the numerical values of nil 
and ni2 for both QD-NH and QD-IH cases, using the observational values of 
| A m| 3 | and Am^. Table-1 gives the calculated numerical values for both 
models [9]. 

In the next step we parameterize the mass matrix in eq.(l) into three 
different types [5,10]: 

Type IA with D=Diag(m!, — m 2 , ni 3 ). The mass matrix of this type can be 
parameterized using two parameters (e,^): 



I e — 2i] —ce —ce 
-ce \-dr] -\- 
\ —ce — \ — dr] 



ce \- dr] -\ - rj \ m 3 . (8) 



This predicts the solar mixing angle, 
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1 llTTAT'ATIT 

umerenc 


QD-NH 


QD-IH 


V\ «~* T" «~* VV1 f\ ■ f'\ T'» O 

pd.1 dllltitt;! a 










c 


1.0 


1.0 


1.0 


1.0 


d 


1.0 


1.0 


1.0 


1.0 


m 3 


0.10 


0.10 


0.08 


0.08 


e 


U. 57972 


0.0015 


0.78004 


0.00169 


77 


0.14602 


0.0649 


0.19628 


-0.08546 


mi (eV) 


0.08674 


0.08675 


0.09340 


0.09340 


m 2 (eV) 


-0.08717 


0.08717 


-0.09380 


0.09380 


m3(el / ) 


0.10 


0.10 


0.08 


0.08 


E \m\eV 


0.27 


0.274 


0.267 


0.274 


Am 2 21 eV 2 


7.6 x lO" 5 


7.6 x 10" 5 


7.6 x 10~ 5 


7.6 x 10~ 5 


\Am 2 23 \eV 2 


2.2 x 10- 3 


2.4 x 10- 3 


2.4 x 10~ 3 


2.4 x 10- 3 


tan 2 ^12 


0.50 


0.50 


0.50 


0.50 


\m ee \ eV 


0.08688 


0.0869 


0.09354 


0.09354 



Table 2: Predictions for tan 2 9i 2 = 0.50 and other parameters consistent with 
observations. 



When we choose the constant parameters c = d — 1.0, we set for the tri- 
bimaximal mixings (TBM) tan20 l2 = -2y/2 (which is tan 2 # 12 = 0.50) and 
the values of e and r\ are calculated for both QD-NH and QD-IH models, by 
using the values of obsrvational data [8] in these two eigenvlue expressions: 
m 1 = (2e— 2r])m 3 , and m 2 = (—e—2r])m 3 obtained for TBM solution in eq.(8). 
The results are given in Table-2 for tan 2 9i 2 = 0.50. The solar mixing angle 
can be further lowered by taking the values c < 1 and d > 1, while retaining 
the earlier values of e and 77 extracted for TBM solution. For tan 2 Q\ 2 = 0.45, 
case the results are shown in Table-3. 

Type-IB with D = Diag (mi,m 2 ,m 3 ): This class of quasi-degenerate mass 
pattern is given by the mass matrix, 



m LL = 



( 1 — e — 2i] ce ce \ 

ce 1 — dr) —i] 
\ ce —i] 1 — drj J 



m 3 . (10) 



This predicts the solar mixing angle, 
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1 llTTAT'ATIT 

umerenc 


QD-NH 


QD-IH 


parameters 


Tvnn TA 


x y \-) u — ± j_> 






c 


0.868 


0.945 


0.868 


0.96 


d 


1.025 


0.998 


1.0 


1.002 


m 3 


0.10 


0.10 


0.08 


0.08 


e 


0.6616 


0.00145 


0.88762 


0.00169 


77 


0.1655 


0.06483 


0.22317 


-0.08546 


mi (eV) 


0.0876 


0.08676 


0.09392 


0.09341 


m 2 (eV) 


-0.0880 


0.08717 


-0.09432 


0.09381 


m 3 (eF) 


0.0996 


0.10002 


0.08 


0.080014 




0.274 


0.274 


0.268 


0.267 




7.7 x 10- 5 


7.3 x 10- 5 


7.6 x 10~ 5 


7.4 x 10~ 5 


|Am2 3 |e\/ 2 


2.2 x 10~ 3 


2.4 x 10~ 3 


2.4 x 10" 3 


2.4 x 10~ 3 


tan 2 9i2 


0.45 


0.45 


0.45 


0.45 


\m ee \ eV 


0.0877 


0.08688 


0.09403 


0.09354 



Table 3: Predictions for tan 2 6*12 = 0.45 and other parameters consistent with 
observations. 

which gives the TBM solar mixing angle with the input values c = 1 and 
d = 1. When e = 0, i] = 0, this leads to rn^ = diag(l, 1, l)m 3 . Like 
in Type-IA, here e and rj values are computed for QD-NH and QD-IH, by 
using observational data [8] in the corresponding two eigenvalue expressions: 
mi = (1 — 2e — 2r])m 3 and m 2 = (1 + e — 2r})m 3 for TBM solution in eq.(10). 

Type-IC with D = Diag(mi,m 2 , — m 3 ): It is not necessary to treat this 
model separately as it is similar to Type-IB except with the interchange of 
two matrix elements (77122) and (77123) m the mass matrix in eq.(10), and this 
effectively imparts an additional odd CP-parity on the third mass eigenvalue 
m 3 in Type-IC. Such change does not alter the predictions of Type-IB. Tables 

2 and 3 present our numerical results for tan 2 9i 2 = 0.50 and tan 2 6i 2 = 0.45 
respectively which are consistent with cosmological bound. 

3 Predictions of baryon asymmetry 

We now apply the above Quasi-degenerate neutrino mass matrices with the 
values of the parameters already fixed in section 2, in the calculation of 
baryon asymmetry of the universe. This gives the second stage for the dis- 
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Type (m,n) M 1 (GeV) M 2 (GeV) M 3 (GeV) 

NH-IA (6^2) 4.8659 x 10 8 -3.5068 x 10 12 9.1256 x 10 14 

(8,4) 3.9414 x 10 6 -3.9774 x 10 10 6.0097 x 10 13 

NH-IB (0) 1.8117 x 10 s 2.6219 x 10 12 3.2528 x 10 14 

(8,4) 1.4994 x 10 6 2.1238 x 10 10 3.2527 x 10 14 

IH-IA 0) 4.5568 x 10 8 -2.8771 x 10 12 1.3273 x 10 14 

(8,4) 3.6910 x 10 6 -2.3241 x 10 10 1.8289 x 10 14 

IH-IB 0) 1.7153 x 10 8 2.8234 x 10 12 3.5098 x 10 14 

(8,4) 41.393 x 10 6 2.4809 x 10 10 3.8091 x 10 14 



Table 4: Heavy right-handed Majorana neutrino mass Mj for QDN with 
normal and inverted ordering mode for tan 2 9\2 = 0.45, using neutrino mass 
matrices given in section 2. The entry (m,n) indicates the type of Dirac 
neutrino mass matrix, as explained in the text. 



crimination among the QDN models. 

For the numerical calculation of baryon asymmetry, we refer to all the 
relevant expressions given in the Appendix. First we translate mass matrices 
rriLL via the inversion of the seesaw formula [11], M RR = —m\ R m~l^ L mi jR - 
We choose a basis for U R where M% R 9 = U^M RR U R = diag(M u M 2 , M 3 ), 
with real and positive mass eigenvalues. We then transform diagonal form 
of Dirac neutrino mass matrix, rriL R = diag(X m , A n , V)v to the U R basis: 
m,L R — > m LR = mL R U R Q where Q = diag(l,e ia ,e l13 ) is the complex matrix 
containing CP-violating Majorana phases introduced by hand. Here A is the 
Wolfenstein parameter and the choice (m,n) in mm gives the type of Dirac 
neutrino mass matrix. The value of the vev is taken as v = 174 GeV. 

At the moment we consider phenomenologically two possible forms of 
Dirac neutrino mass matrices such as (i) (m,n) = (6,2) for the charged- 
lepton type mass matrix, and (ii) (8,4) for up-quark type mass matrix. In 

this prime basis the Dirac neutrino Yukuwa coupling becomes h = ^f - which 
enters in the expression of CP-asymmetry e in (14) and (23) given in the 
Appendix. The new Yukawa coupling matrix h also becomes a complex 
quantity, and hence the term Im(h^h)ij appearing in lepton asymmetry e±, 
gives a non-zero contribution. A straight forward simplification shows that 
(h)h)\j = (Q22) 2 Q22^2 + (Ql 1 ) 2 Ql 3 R2 where R 2 ,3 are real parameters. After 
inserting the values of phases, the above expression leads to Im{h)h)\j = 
— [R 2 sin2(a — (3) + R 3 sin2a] which imparts a non-zero CP asymmetry for 
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lype 


(m,n) 






e l 




VlB 




V3B 




NH-IA 


(6,2) 


3.73 x 10" 


-6 


1.92 x 10- 


-7 


9.07 x 10- 


-12 


2.11 x 10- 


-li 




(8,4) 


3.03 x 10- 


-8 


1.55 x 10" 


-9 


7.32 x 10" 


-14 


1.71 x 10" 


-13 


NH-IB 


(6,2) 


5.31 x 10" 


-7 


1.12 x 10- 


-14 


1.41 x 10- 


-18 


5.67 x 10" 


-13 




(8,4) 


4.30 x 10" 


-9 


8.87 x 10" 


-17 


1.12 x 10" 


-20 


4.71 x 10" 


-15 


IH-IA 


(6,2) 


3.77 x 10" 


-6 


1.94 x 10" 


-7 


8.50 x 10" 


-12 


1.98 x 10- 


-11 




(8,4) 


3.05 x 10" 


-8 


1.57 x 10" 


-9 


6.88 x 10" 


-14 


1.60 x 10- 


-13 


IH-IB 


(6,2) 


5.31 x 10" 


-7 


9.75 x 10" 


-15 


1.15 x 10" 


-18 


5.95 x 10" 


-13 




(8,4) 


4.30 x 10" 


-9 


7.80 x 10" 


-17 


9.17 x 10" 


-21 


4.76 x 10" 


-15 



Table 5: Values of CP asymmetry e\ and baryon asymmetry [j]\BiVzb) for 
all quasi-degerate models, with tan 2 Qyi = 0.45, using neutrino mass matrices 
given in the text. The entry (m,n) in indicates the type of Dirac neutrino 
mass matrix taken as charged lepton mass matrics (6,2) or up quark mass 
matrix (8,4) as explain in the text. 



particular choice of (a, (3). 

In our numerical estimation of lepton asymmetry, we chose some arbitrary 
values of a and (3 other than | and 0. For example, light neutrino masses 
(mi, -m 2 , m 3 ) of Type-IA model, leads to M%%? = diag(M 1 , -M 2 , M 3 ), and 
we thus fix the Majorana phase Q = diag(l, e ia , e t/3 ) = diag(l, e l ( 7r / 2 + 7I 7 4 ) ) e l7r//4 ) 
for a = (7r/2 + 7r/4) and (3 = tt/4. The extra phase 7r/2 in a absorb the neg- 
ative sign before heavy Majorana mass. In our search programmes such a 
choice of the phase leads to highest numerical estimation of lepton CP asym- 
metry. 

In Table 4 we give numerical prediction on three heavy right-handed 
Majorana neutrino masses from these neutrino mass models under consider- 
ation for the case of tan 2 9 12 = 0.45. The three heavy right-handed Majorana 
mass matrices which are constructed through the inversion of Type-I seesaw 
mechanism, for two choices of diagonal Dirac-neutrio mass matrix discussed 
before. The corresponding baryon asymmetries t]b are estimated for both 
non-flavour t]ib and flavour tj^b leptogenesis respectively in Table 5. As ex- 
pected, there is an enhancement in baryon asymmerty by a magnitude of 
order one (approximately) when flavour dynamics is included (see Table 5 
in Type-IA model with charged lepton mass matrix) [12,13]. Type-IA with 
charged lepton Dirac neutrino mass matrix, is the only model sensitive and 



8 



Type (m,n) T% m < T H < T% ax (GeV) M^ m < < M^ ax (GeV) 

NH-IA (6,2) 5.51 x 10 6 < T R < 4.87 x 10 7 9.73 x 10 9 < < 8.59 x 10 10 
(8,4) 2.21 x 10 2 < T R < 3.94 x 10 4 7.88 x 10 6 < M$ < 1.40 x 10 9 

NH-IB (0) 3.7 x 10 8 < T R < 1.85 x 10 6 3.70 x 10 8 < M$ < 1.89 x 10 2 
(8,4) 2.96 x 10 13 < T R < 1.49 x 10 4 2.99 x 10 6 < < 1.52 x 10~ 2 

IH-IA {0} 5.12 x 10 5 < T R < 4.56 x 10 6 9.11 x 10 8 < M$ < 8.11 x 10 9 
(8,4) 5.12 x 10 5 < T R < 3.69 x 10 4 7.38 x 10 6 < < 5.32 x 10 5 

IH-IB (6,2) 3.84 x 10 12 < T R < 1.72 x 10 7 3.44 x 10 8 < < 1.54 x 10 2 
(8,4) 3.88 x 10 12 < T R < 1.39 x 10 4 2.79 x 10 6 <M <S> < 9.99 x 1(T 3 

Table 6: Theoretical bound on reheating temperature Tr and inflaton masses 
Mff, in non-thermal leptogenesis, calculated using data from Table 5, for all 
neutrino mass models with tan 2 9\2 = 0.45. 

consistent with data on observed baryon asymmetry. 

In case of non-thermal leptogenesis, the lightest right-handed Majorana 
neutrino mass Mi from Table 4 and the CP asymmetry t\ from Table 5, 
for all the neutrino mass models, are used in the calculation of the bounds: 
T min < t r < T^ ax and M% in < < M^ ax in Table 6. The baryon 
asymmetry Y B — ^ = 8.7 x 10~ n is taken as input value from WMAP 
obseervational data. Only those neutrino mass models which simultaneously 
satisfy the two constraints, T^ ax > T r R nin and M^ ax > M^ in , could survive in 
the non-thermal leptogenesis scenario. The predicted inflaton mass as ~ 
10 11 GeV for reheating temperature Tr = 10 6 GeV are needed to produce 
the observed baryon asymmetry of the universe [14,15]. From Table 6, the 
neutrino mass nodels with (m,n) which are compatible with ~ (10 10 — 
10 13 ) GeV and T R w (10 6 -10 7 ) GeV are listed as NH-IA (6,2) and IH-IA (6,2) 
only. Again in order to avoid gravitino problem [8] in supersymmetric models, 
one has the bound on reheating temperature, Tr (10 6 — 10 7 ) GeV. This 
implies that Type-IA where Dirac neutrino mass matrix is taken as charged 
lepton mass matrix is the only model consistent with the observed baryon 
asymmetry. These findings nearly agree with flavoured thermal leptogenesis 
for Type NH-IA (6,2) model in Table 5. This result is also consistent with 
quasi-degenerate in inverted hierachical (IA) models for charged lepton mass 
matrix (6,2) [15]. 
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4 Summary and Conclusions 



We have studied the effects of Majorana CP phases on the prediction of ab- 
solute neutrino masses in three types of QDN models having both normal 
and inverted hierarchical patterns within /i-t symmetry. These predictions 
are consistent with data on the mass squared difference derived from vari- 
ous oscillation experiments, and from the upper bound on absolute neutrino 
masses in Ov/3/3 decay as well as cosmological upper bound of J2i m i < 0.28 
eV. It has been found that QDN models with mi < 0.09 eV, are still far from 
discrimination and hence the quasi-degenerate pattern is not yet ruled out. 
The prediction on solar mixing angle is also found to be lower than TBM 
value viz, tan 2 12 =0.45 which coincides with the best-fit in the neutrino 
oscillation data [8]. 

In the next stage, left-handed Majorana neutrino mass matrices tjill have 
been employed to estimate the baryon asymmetry of the universe, in both 
thermal and non-thermal leptogenesis scenario (Tables 5-6). We use the CP 
violating Majorana phases derived from right-handed Majorana mass matrix, 
and also Dirac neutrino mass matrix as either charged lepton mass matrix or 
up-quark mass matrix. We also observe some enhancement effects in flavour 
leptogenesis [20] compared to non-flavour leptogenesis by a magnitude of or- 
der one. The predicted inflaton mass needed to produce the observed baryon 
asymmetry of the universe is found to be ~ 10 10 GeV corresponding to 
the reheating temperature Tr = 10 6 GeV [14,15]. The analysis shows that 
quasi-degenerate model in normal hierarchical pattern (NH-IA) with charged 
lepton Dirac mass matrix (6,2), appears to be a favourable choice in nature. 
The quasi-degenerate model in inverted heirarchy(IH-IA) with the charged 
lepton Dirac mass matrix (6,2), is not completely ruled out. The results pre- 
sented in this article are new and have subtle hints in the discrimination of 
neutrino mass models. This could establish the quasi-degenerate neutrinos 
as natural physical neutrinos in the neutrino oscillation experiments. 
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Appendix 



Thermal leptogenesis 

The light left-handed Majoarana neutrino mass matrix ttill and heavy right- 
handed Majorana mass matrix M RR are related through the canonical seesaw 
formula (known as Type I) [11] in a simple way: 

m LL = -m LR M RR m T LR (12) 

where m^ R is the Dirac neutrino mass matrix. For our calculation of lepton 
asymmetry, we consider the model [16,17] where the asymmetry decay of 
the lighest of the heavy right-handed Majorana neutrinos, is assumed. The 
physical Majorana neutrino N R decays into two modes: N R — > Il + <P, N r — > 
Il + ip where II is the lepton and II is the antilepton; (p and (p are the Higgs 
and anti-Higgs particles respectively. The braching ratio for these two decay 
modes are likely to be different. The CP asymmetry which is caused by the 
interference of tree level with one-loop corrections for the decays of the ligh 
test of heavy right-handed Majorana neutrino N R is defined by [16,18] 

r - f 

e = 7^ (13) 

r + r y J 

Here T = T(Ni — > l L (p) and T = T(Ni — > l L (p) are the decay rates. A 
perturbative calculation from the interference between tree level and vertex 
plus self energy diagrams, gives [19] the lepton asymmetry e\ for non-SUSY 

case as 

« = ^(4, s/ mlih,h)M U u U )] (14) 

where f(x) and g(x) represent the contributions from vertex and self energy 
corrections respectively, f(x) = \/x[— 1 + (x + l)/n(l + ^)],g(x) = For 
hierarchical right-handed neutrino masses where x is large, we have the ap- 
proximation [16], f(x) + g(x) = This simplifies to 

_ 3 Im^hf^ M, Im[(h^h)l 3 ] M 1 
1 16tt [ (tfh)u M 2 (h!h)u M 3 J 1 ' 

where h = is the Yukawa coupling of the Dirac neutrino mass matrix in 
the diagonal basis of M RR and v= 174 GeV is the vev of the standard model. 
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In term of light Majorana neutrino mass matrix m LL the above expression 
can be simplified to 

ei = -li^M(^W^*)n] (16) 

For quasi-degenerate spectrum of the heavy right-handed Majorana neu- 
trino masses, i.e., for Mi ps M 2 < M 3 the asymmetry is largely enhanced by 
a resonance factor and in such situation, the lepton asymmetry is modified 
[21] to 

1 Im[{h)h)\ 2 ] 
£l = 8^ {Wh) u R (17) 

where R = ^'^}^2+r^M' 1 anc ^ ^ 2 = ^ h ^ M2 . It can be noted that in case 
of SUSY, the functions f(x) and g(x) are given by f(x) = y/xln^l + -) and 

9( x ) = fzfi f° r large x one can have f(x) + g(x) ps 3/y/x. Therefore the 
factor 3/8 will appear in place of 3/16 in the expression of CP asymmetry 
[20] in eq. (15). The CP asymmetry parameter 6j is related to leptonic 
asymmetric parameter Y L as 

y L « HI^L = j- e -^l (18) 

where ni is the lepton number density, ul is the anti- lepton number density, 
s is the entropy density, Ki is the dilution factor for the CP asymmetry e^, 
and g*i is the effective number of degrees of freedom at temperature T — Mi. 
The baryon asymmetry Yb produced through the sphaleron transition of Yl 
while the quantum number B-L remains conserved, is given by [22] 

F B = — = CYb-l = CY L (19) 
s 

where C = 2 2Nf+i3N h ' is the number of fermion families and Nh is the 
number of Higgs doublets. Since s = 7.04n 7 the baryon number density over 
photon number density n 7 corresponds to the observed baryon asymmetry 
of the universe [23], 

= (— ) SM « (20) 

where d ^ 0.98 x 10 -2 is used in the present calculation. In case of MSSM, 
there is no major numerical change with respect to the non-supersymmetric 
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case in the estimation of baryon asymmetry. One expects approximate en- 
hancement factor of about \/2y/2 for strong (weak) washout regime [20]. 

In the expression for baryon-to-photon ratio in eq. (20), K\ describes 
the washout factor of the lepton asymmetry due to various lepton number 
violating processes. This efficiency factor (also known as dilution factor) 
mainly depends on the effective neutrino mass 

_ (h*h) u u 2 

mi = (21) 
where v is the electroweak vev, v = 174 GeV. For lO~ 2 eV <rh\< 10 3 e^, 
the washout factor K\ can be well approximately by [19,24] 

Ki(ii)=a3 C |N ^ r °' 6 (22) 

The value of K\ is valid only for the given range of m [25]. 
Flavoured thermal leptogenesis 

It is inevitable to include the flavour effects in thermal leptogenesis [26] 
and study its effects on the enhancement in baryon asymmetry over the 
single flavour approximation. In the flavour basis the equation for lepton 
asymmetry in N\ — > l a ip decay where a = (e, //, r) becomes 

€aa = JT (hU,\ t ^ ^[^i(^ t ^)ij^i]^(^)+E /m [^i(^ t ^)jiM n —- d 
8tt {Wh)u i=2 3 . (1 - Xj) 

(23) 

of-equilibrium situation is given by K a = S^-. Here 8n ^2 ~ 1-1 x 10" 3 eV, 
and m aa = ^V 2 . This leads to the baryon asymmetry of the universe, 

mB = 5£ „ io- 2 £ e aa K a ~ 10- V £ ^ (24) 

a a m aa 

For single flavour case, the second term in e aa vanishes when summed over 
all flavours. Thus 

I i 

d = E e «« = (XtX) ? M(^)ijM^)- ( 25 ) 
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Here we have x 3 - = ^4 and <7 (£_,-) ~ |~^- The efficiency factor for the out- 



This leads to baryon asymmetry, 



7] 1B « 10 _2 m*^ = 10" 2 /tiei (26) 
m 



where e 1 = £ Q e aa and m = £ a m c 



Non-themal leptogenesis 

We now focus our attention on the non-thermal leptogenesis scenario [27] 
where the right-handed neutrinos are produced through the direct non-thermal 
decay of the inflaton 0. We follow the standard procedure [28] where non- 
thermal leptogenesis and baryon asymmetry in the universe had been studied 
in different neutrino mass models whereby some mass models were excluded 
using bound from below and from above on inflaton mass and reheating 
temperature after inflation. We start with the inflaton decay rate given by 

r, = r|HW»g (27) 

where Aj are the Yukawa coupling constants for the interaction of three heavy 
right-handed neutrinos iVj with the inflaton <p of mass M^. The reheating 
temperature after inflation is given by the expression, 

T * = (^)*F* M p) h (28) 

where M p 2.4 x 10 18 Gev is the reduced Plank mass [29] and g* is the 
effective number of relativistic degrees of freedom at reheating temperature. 
For SM we have g*= 106.75, and for MSSM g m = 228.75. If the inflaton 
dominantly couples to N iy the branching ratio of this dacay process is taken 
as BR ~ 1, and the produced baryon asymmetry of the universe can be 
calculated by the following relation [30], 

Ys = ? = CY L = (29) 

where Y L is the lepton asymmetry genarated by CP-violating out-of- 
equilibrium decays of heavy neutrino Ni and Tr is the reheating temper- 
ature. The fraction C has the value C=-28/79 for SM and C=-8/15 in the 
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MSSM. The observed baryon asymmetry measured in WMAP data, r\ B = 
2a = 6.5 x KT 10 [3], where s=7.04 n 7 is ralated to Y B = n B /s = 8.6 x KT 11 
in eq.(29). From eq.(29) the connection between T B and is expressed as, 

Tr = (^-W, (30) 

The above expression is supplemented by two more boundary conditions [28] : 
(i) lower bound on inflaton mass, > 2Mi coming from allowed kinemat- 
ics of inflaton decay to two right-handed Majorana neutrinos Ni, and (ii) 
an upper bound for the reheating temperature, Tr < O.OlMi coming from 
out-of-thermal equilibrium decay of JVi. Using the observed central value of 
the baryon asymmetry Y B and theoretical prediction of CP asymmetry ei 
in eq.(30), one can establish the relation between T R and for each neu- 
trino mass model. The lightest right-handed neutrino mass Mi and the CP 
asymmetry e\ neutrino mass models are used in the calculation of theoretical 
bounds: T^ in < T R < T^ ax and M^ in < < M^ ax following eq.(19) 
along with other two boundary conditions cited above. Only those models 
which satisfy the constraints T R ax > T R m and M™ n < M™ ax simultneously 
can survive in the non-thermal leptogenesis. 
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